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The  strain  held  of  isolated  screw  and  edge  dislocation  cores  in  aluminum  are  calculated  using  density- 
functional  theory  and  a  flexible  boundary  condition  method.  Nye  tensor  density  contours  and  differential 
displacement  fields  are  used  to  accurately  bound  Shockley  partial  separation  distances.  Our  results  of  5- 
7.5  A  (screw)  and  7. 0-9. 5  A  (edge)  eliminate  uncertainties  resulting  from  the  wide  range  of  previous 
results  based  on  Peierls-Nabarro  and  atomistic  methods.  Favorable  agreement  of  the  predicted  cores  with 


limited  experimental  measurements  demonstrates 

location  cores. 
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Plastic  deformation  in  metals  is  controlled  by  the  crea¬ 
tion,  motion,  and  evolution  of  line  defects  known  as  dis¬ 
locations.  Dislocations  move  at  stresses  far  below  that 
required  to  irreversibly  deform  a  defect-free  crystal  [1]. 
They  exhibit  multiscale  behavior  characterized  by  an 
atomic-length  scale  “core”  at  the  dislocation  center  and 
a  long-range  strain  held  that  decreases  as  the  inverse  of  the 
distance  from  the  core.  In  fee  metals  dislocations  can 
reduce  elastic  energy  by  separating  into  Shockley  partial 
dislocations  connected  by  a  stacking  fault  [1].  For  materi¬ 
als  with  large  stacking  fault  energies  (i.e.,  aluminum)  this 
separation  can  be  quite  narrow  making  experimental  in¬ 
vestigation  of  the  core  structure  difficult  [2].  Several  first- 
principles  methods  have  been  developed  that  attempt  to 
address  the  coupling  of  long-range  elastic  distortion  to  the 
local  displacements  in  the  core  region  [3].  These  and 
atomistic  calculations  have  led  to  a  wide  range  in  the 
predicted  structure  of  dislocations  in  A1  (cf.  Table  II). 

Accurate  modeling  of  the  dislocation  core — the  local 
strain  held  produced  by  the  partials  and  their  separation 
distance — is  key  to  predictive  modeling  of  plastic  behav¬ 
ior  in  metals.  For  example,  the  motion  of  dislocations  in 
pure  A1  under  stress  is  governed  by  geometry-dependent 
interactions  with  sound  waves  (phonon  drag)  and  multiple 
slip-planes  (cross-slip)  [1].  Recent  models  of  dynamic 
strain  aging  in  Al-Mg  alloys  show  that  solute  binding 
energies  and  diffusion  activation  enthalpies  within  the 
dislocation  core  region  control  related  macroscopic  phe¬ 
nomena  [12,13].  The  most  familiar  of  these  is  the  Porte vin- 
LeChatelier  (PLC)  effect  which  is  commonly  found  in 
commercial  Al-Mg  alloys  [14].  Accurate  quantitative  in¬ 
puts  to  these  theories  are  unavailable  due  to  the  approx¬ 
imations  inherent  in  atomistic  methods,  and  the  boundary 
conditions  used  in  previous  first-principles  approaches.  For 
many  alloys  correctly  describing  the  atomic  scale  pro¬ 
cesses  controlling  macroscopic  phenomena  will  require 
quantum  mechanical  predictions  with  realistic  boundary 
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conditions,  this  is  particularly  true  for  the  A1  alloys.  These 
alloys  have  a  complex  chemistry  and  flow  behavior  and 
play  a  vital  role  in  a  myriad  of  technical  applications. 

In  this  Letter  we  employ  electronic  structure  calcula¬ 
tions  based  upon  density-functional  theory  (DFT)  and  a 
flexible  boundary  condition  method  to  optimize  the  strain 
field  of  isolated  screw  and  edge  dislocations  in  pure  Al. 
While  this  methodology  has  been  applied  to  dislocations  in 
transition  metals  such  as  bcc  molybdenum  and  tantalum 
[15],  fee  iridium  [16],  and  ordered  alloys  [17],  the  present 
treatment  is  the  first  application  to  dislocations  in  a  free 
electron  metal.  To  treat  the  screw-  and  edge-character 
dislocation  geometries  accurately,  simulation  cells  are  de¬ 
signed  to  minimize  electrostatic  and  strain  field  interac¬ 
tions  with  the  required  surface  boundary.  Moreover,  the 
long-range  dislocation  strain  field  is  treated  self- 
consistently  using  the  first-principles  lattice  Green  func¬ 
tion  boundary  condition  method  (FP-GFBC)  [15].  This 
DFT  technique  correctly  couples  the  strain  field  of  the 
dislocation  core  to  the  long-range  elastic  field.  When  ap¬ 
plied  to  Al  the  method  produces  a  separation  distance  for 
the  edge-character  dislocation  which  is  in  good  agreement 
with  the  available  experimental  measurements  of  these 
dislocation  cores  [4]. 

Computational  methods .  — The  ab  initio  calculations  are 
performed  with  VASP  [18,19],  a  density-functional  code 
using  pseudopotentials  and  a  plane- wave  basis.  For  Al 
we  treat  the  3  s  and  3  p  states  as  valence  and  employ  a 
plane- wave  kinetic-energy  cutoff  of  192  eV  to  ensure  an 
accurate  treatment  of  the  valence  electrons.  An  initial  study 
was  undertaken  to  determine  the  optimum  pseudopotential 
representation.  Structural  parameters  are  calculated  at 
three  levels  of  increasing  accuracy:  ultrasoft  pseudo¬ 
potentials  [20,21]  (USP)  in  the  local  density  approxima¬ 
tion,  USP’s  in  the  generalized  gradient  approximation 
(GGA)  of  Perdew  and  Wang  [22]  and  finally  projected 
augmented  wave  (PAW)  method  within  GGA.  The 
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&-point  meshes  for  these  initial  calculations  are  15  X  15  X 
15,  while  the  dislocation  calculations  employ  a  1X1X8 
mesh,  both  use  a  Methfessel-Paxton  smearing  of  1  eV. 
Table  I  shows  the  calculated  parameters  including  the  shear 
modulus  [G  =  (C44  +  Cn  —  C12)/3]  relating  stress  and 
strain  along  (110)  in  the  (111)  plane.  The  lattice  constant, 
C44  and  G  significantly  improve  in  going  to  US-GGA  and 
then  PAW-GGA.  The  wide  range  of  measured  stacking 
fault  energies  reflects  the  difficulty  of  resolving  Shockley 
partial  dislocations  in  Al.  The  current  results  are  smaller 
than  that  calculated  previously  using  norm-conserving 
pseudopotentials  (161  mJ/m2),  and  better  match  aug¬ 
mented  plane-wave  results  (126  mJ/m2)  [27,28].  As  a 
compromise  between  efficiency  and  precision  the  US- 
GGA  potentials  are  used  in  the  dislocation  core 
calculations. 

Calculations  of  an  isolated  dislocation  requires  careful 
construction  of  the  simulation  cell.  The  cut  in  the  lattice 
produced  by  the  dislocation  (screw  or  edge)  is  incompat¬ 
ible  with  the  periodic  symmetry  imposed  by  the  supercell, 
so  in  these  calculations  a  free  surface  is  introduced  to 
contain  the  dislocation.  The  vacuum  region  produces  a 
charge  dipole  at  the  surface  with  slowly  decaying  Friedel 
oscillations  in  the  charge  density.  Our  reference  calcula¬ 
tions  show  that  the  atomic  forces  resulting  from  the  charge 
dipole  extend  up  to  10  A  from  the  surface.  The  relaxation 
of  the  edge-character  dislocation  is  calculated  in  a  58  A  X 
58  A  X  4.945  A  periodic  cell,  using  552  atoms  in  a  24.2  A 
radius  cylindrical  slab  surrounded  by  vacuum.  The  inner  38 
atoms  (region  I)  around  the  core  are  relaxed  using  ab  initio 
forces,  while  the  ab  initio  forces  of  the  next  99  atoms 
(region  II)  are  relaxed  using  the  lattice  Green  function. 
The  positions  of  the  outer  415  atoms  (region  III)  are 
updated  when  minimizing  the  forces  in  region  II  [15]. 
The  relaxation  of  the  screw-character  dislocation  is  calcu¬ 
lated  in  a  68.6  A  X  68.6  A  X  2.857  A  periodic  cell,  using 
437  atoms  in  a  28.3  A  radius  cylindrical  slab  surrounded  by 
vacuum.  The  inner  63  atoms  are  in  region  I,  the  next  66 
atoms  in  region  II,  and  the  remaining  308  atoms  in 
region  III.  For  both  dislocations,  the  ab  initio  forces  are 
relaxed  iteratively  using  conjugate-gradient  minimization 
in  region  I,  followed  by  lattice  Green  function  updates  to 
minimize  the  forces  in  region  II,  until  the  atomic  forces  are 


TABLE  I.  Calculated  and  measured  lattice  constant,  elastic 
constants  (Mbar),  and  stacking  fault  energy  (mJ/m2)  for  fee 
aluminum  [23-26].  The  VASP  calculated  values  are  based  on 
different  levels  of  approximation  of  the  pseudopotentials  and  the 
local  density  approximation  as  described  in  the  text. 


Method 

a  [A] 

C„ 

C12 

C44 

G 

7 

US-LDA 

3.9780 

1.218 

0.634 

0.373 

0.319 

134 

US-GGA 

4.0413 

1.085 

0.568 

0.305 

0.274 

124 

PAW-GGA 

4.0488 

1.0593 

0.5708 

0.3078 

0.264 

122 

Experiment 

4049 

1.130 

0.666 

0.278 

0.247 

120-200 

less  than  5  meV  A-1.  Convergence  of  the  strain  field  with 
cell  size  was  verified  by  optimizing  larger  simulation  cells. 

A  lattice  Green  function  for  Al  is  derived  using  the 
calculated  elastic  constants  and  the  force-constant  matrix 
[29].  The  force-constant  matrix  is  calculated  using  the 
direct  force  method  in  a  3  X  3  X  3  simple  cubic  supercell 
(108  atoms),  with  a  8  X  8  X  8  &-point  mesh.  The  estimated 
relative  error  [29]  of  the  calculated  lattice  Green  function 
due  to  finite  size  errors  is  1.2%. 

Results .  — Figure  1  shows  the  differential  displacements 
(DD)  and  Nye  tensor  distribution  for  the  screw  and  edge 
components  of  the  a/ 2(110)  screw-character  dislocation 
including  measurements  of  partial  separation.  Differential 
displacements  show  the  relative  displacement  from  a  per¬ 
fect  lattice  of  two  nearest  neighbor  atoms  projected  in  the 
screw  or  edge  direction  by  a  vector  of  the  corresponding 
length  between  the  two  atoms.  The  partial  dislocations  are 
identified  by  a  triad  of  atoms  whose  three  vectors  have  a  net 
chirality  (either  clockwise  or  counterclockwise).  The  Nye 
tensor  distribution  uses  the  atomic  displacements  of  an 
atom  relative  to  each  nearest  neighbor  to  produce  a  con¬ 
tinuous  distribution  of  Burgers  vector  components;  the 
individual  values  are  projected  onto  the  plane  normal  to 
the  dislocation  line  and  smoothly  interpolated  [30].  In  this 
case  the  partial  dislocations  are  identified  by  local  extrema 
in  the  Nye  tensor  distribution.  The  width  of  the  Nye  tensor 
distribution  shows  the  partial  dislocations  spread  over 
several  lattice  positions,  which  overlap  for  the  screw  com- 
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FIG.  1  (color).  Atomic  core  structure  (screw  and  edge  compo¬ 
nents)  of  the  two  partials  of  the  a/ 2[1 10]  screw-character  dis¬ 
location.  The  atomic  positions  (circles)  are  superimposed  on  the 
differential  displacements  (arrows)  and  the  Nye  tensor  distribu¬ 
tion  (color  map);  see  text  for  details.  Using  the  screw  compo¬ 
nents,  the  triad  of  arrows  with  a  net  chirality  identifies  the 
position  of  the  partials  whose  separation  is  d®D.  From  the 
edge  components,  the  centers  of  the  partials  are  identified  as 
the  extrema  in  the  Nye  tensor  distribution  whose  separation  is 
d^ye.  Using  these  two  measurements  the  partial  core  separation 
is  shown  to  be  in  the  range  of  5.0 -7.5  A. 
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ponent.  The  DD  plot  of  the  screw  component  gives  clear 
indications  of  the  partial  dislocations,  as  do  the  extrema  in 
the  Nye  tensor  distribution  of  the  edge  components.  Using 
these  two  approximations  the  partial  core  separation  is 
shown  to  range  from  d®D  =  7.5  A  to  <ifye  =  5.0  A. 
Comparing  the  screw  and  edge  components  we  see  that 
separation  distances  predicted  by  each  method  are  self- 
consistent;  however,  the  Nye  tensor  gives  slightly  narrower 
separation  distance.  The  range  of  values  from  the  same 
geometry  illustrates  the  difficulties  of  defining  the  center  of 
the  partial  dislocations  to  a  resolution  better  than  the  lattice 
spacing. 

Figure  2  shows  the  edge  and  screw  components  of  the 
edge-character  dislocation  which  has  a  broader  partial 
separation  compared  with  the  screw-character  dislocation 
(Fig.  1).  The  separation  of  partials  is  computed  in  an 
identical  fashion  as  for  the  screw-character  dislocation 
and  in  this  case  the  DD  predicts  d =  9.5  A  and  the 
Nye  tensor  distribution  predicts  d^yQ  =  7.0  A.  The  nar¬ 
rower  width  of  each  partial,  along  with  the  smaller  lattice 
spacing  in  the  separation  direction  increases  the  accuracy 
of  the  separation  prediction  compared  with  the  screw- 
character  dislocation.  However,  as  with  the  screw  disloca¬ 
tion  the  Nye  tensor  produces  consistently  smaller  partial 
splitting  for  the  edge  dislocation  as  compared  to  the  dif¬ 
ferential  displacement  method.  Accurate  prediction  of  the 
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FIG.  2  (color).  Atomic  core  structure  (edge  and  screw  compo¬ 
nents)  of  the  two  partials  of  the  a/ 2[110]  edge-character  dis¬ 
location.  As  in  Fig.  1  the  atomic  positions  (circles)  are 
superimposed  on  the  differential  displacements  (arrows)  and 
the  Nye  tensor  distribution  (color  map).  Using  the  edge  compo¬ 
nents,  the  triad  of  arrows  with  a  net  chirality  identifies  the 
position  of  the  partials  whose  separation  is  d^D.  From  the  screw 
components,  the  center  of  the  partials  is  identified  as  the  extrema 
in  the  Nye  tensor  distribution  whose  separation  is  d^ye.  Using 
these  two  measurements  the  partial  core  separation  is  shown  to 
be  in  the  range  of  7.0-9.5  A. 


dislocation  geometry  is  a  crucial  first  step  to  understanding 
the  effects  of  stress,  temperature,  and  solutes  on  plasticity 
in  Al. 

Table  II  reviews  observed  and  calculated  separation 
distances  for  the  Al  screw-  and  edge-character  partial 
dislocations,  revealing  a  broad  range  of  results  [4-11]. 
Our  result  for  the  partial  separation  distance  in  the  edge- 
character  dislocation  is  in  good  agreement  with  available 
weak-beam  transmission  electron  microscopy  [4] .  An  ear¬ 
lier  density-functional  theory  calculation  employed  the 
quasicontinuum  DFT  method,  where  a  small  first- 
principles  cell  containing  the  dislocation  is  embedded  in 
the  strain  field  produced  by  atoms  interacting  through  an 
interatomic  force  law  based  on  the  embedded  atom  method 
(EAM)  [5].  The  disagreement  between  the  QC-DFT-EAM 
calculation  and  both  our  calculation  and  the  experimental 
results  most  likely  stems  from  the  small  (84  atom)  cell  used 
for  the  first-principles  calculations.  Our  reference  calcula¬ 
tions  indicate  that  these  cell  dimensions  are  too  small  to 
isolate  the  dislocation  core  elastically  or  electrostatically 
from  the  surface  of  the  metal  cluster.  Peierls-Nabarro 
models  [6,7]  are  not  well  suited  for  describing  details  of 
the  strain  field  in  the  partial  core  as  it  is  represented  only 
indirectly  through  the  generalized  stacking  fault  energy. 
Different  implementations  of  the  Peierls-Nabarro  model 
predict  different  separations,  though  both  are  based  on 
first-principles  evaluations  of  the  generalized  stacking  fault 
energy.  The  results  of  Schoeck  show  good  agreement  with 
the  measured  separation  distance  between  Shockley  parti¬ 
als  for  an  edge  dislocation  [7].  Several  atomistic  calcula¬ 
tions  for  Al  dislocations  are  based  on  central  force 
potentials  [11,31]  and  a  new  angular  dependent  potential 
[10].  Surprisingly,  groups  using  the  same  atomistic  poten- 

TABLE  II.  Shockley  partial  splitting  distances  for  <2/2(110) 
screw  and  edge  dislocations.  Experimental  measurements  for  the 
edge  dislocation  are  from  weak-beam  transmission  electron 
microscopy  (WB-TEM)  [4].  FP-GFBC  is  the  current  work 
showing  two  measurements  of  the  separation  (cf.  Figs.  1  and 
2).  The  quasicontinuum  calculation  (QC-DFT-EAM),  the 
Peierls-Nabarro  methods,  and  the  embedded  atom  methods 
(EAM,  EAM2,  BAM)  are  described  briefly  in  the  text. 


Method 

Screw  [A] 

Edge  [A] 

Experiment 

WB-TEM  [4] 

Density-functional  theory 

8.0 

FP-GFBC(this  work) 

7.5  (DD),  5.0  (Nye) 

9.5  (DD),  7.0  (Nye) 

QC-DFT-EAM  [5] 

Peierls-Nabarro 

5.6 

Semidiscrete  [6] 

2.1 

3.5 

Generalized  [7] 

Atomistic 

7.8 

EAM  [8] 

4.9 

EAM  [9] 

8.8 

14.7 

EAM  [10] 

16.0 

16.0 

EAM2  [11] 

5.2 -5.5  (Nye) 

5.4-6.2  (Nye) 

BAM  [10] 

5.0/15.0 
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tial  [31]  have  reported  dramatically  different  results  for  the 
separation  distance  of  the  screw  dislocation,  ranging  from 
4.9  to  16  A  [8-10];  this  variation  is  possible  the  result  of 
different  treatments  for  the  boundary  conditions.  More 
recently  developed  EAM  potentials  for  A1  from  Mishin 
et  al.  produce  smaller  separation  distances  than  previous 
EAM  potentials,  in  better  agreement  with  the  current  work 

[11].  Finally,  while  elasticity  theory  suggests  only  one 
equilibrium  configuration  for  Shockley  partials  spreading 
on  the  glide  plane,  Srinivasan  et  aV s  potential  produced 
two  stable  partial  separations  at  5.0  and  15  A.  Such  results 
are  consistent  with  strong  short-range  angular  bonding 
terms  in  the  atomistic  potential.  This  broad  range  in  the 
predicted  partial  separation  illustrates  the  sensitivity  of  the 
calculation  to  both  the  treatment  of  the  local  atomic  inter¬ 
actions  and  the  long-range  boundary  conditions.  This  re¬ 
inforces  the  need  for  a  self-consistent  predictive  method 
for  modeling  dislocation  core  structure. 

Conclusion. — We  have  predicted  the  strain  fields  of 
isolated  Al  edge  and  screw  dislocation  using  density- 
functional  theory  combined  with  a  flexible  boundary  con¬ 
dition  method.  The  calculated  Shockley  partial  splitting 
distance  for  the  edge  dislocation  is  in  good  agreement  with 
available  experimental  data.  These  predictions  signifi¬ 
cantly  narrow  the  range  of  splitting  distances  for  both 
screw  and  edge  dislocations  thereby  eliminating  confusion 
stemming  from  the  wide  range  of  Al  values  in  the  litera¬ 
ture.  The  accurate  prediction  of  the  atomic-scale  core 
geometries  for  Al  screw-  and  edge-character  dislocations 
is  needed  to  inform  higher-length  scale  models  of  plasticity 
in  Al  alloys,  as  well  as  other  fee  metals.  Nowhere  is  this 
more  important  than  in  models  of  diffusion  processes  in  Al 
alloys  (e.g.,  Al-Mg)  where  thermodynamic  inputs  must  be 
obtained  from  accurate  solute  or  dislocation  energetics. 
The  wide  range  of  dislocation  geometries  produced  by 
atomistic  simulations  implies  that  such  methods  cannot 
provide  accurate  solute-dislocation  interaction  energies. 
FP-GFBC  optimized  screw  and  core  structures  can  now 
serve  as  the  basis  for  computing  quantities  such  as  the  Al 
core  activation  enthalpy  for  diffusion  and  solute- 
dislocation  energetics  for  such  important  alloying  agents 
as  Mg,  Cu,  Mn,  Li,  Zn,  and  Si,  and  provide  a  basis  for 
understanding  the  influence  of  trace  elements  such  as  Fe. 
While  these  calculations  are  challenging,  precise  predic¬ 
tions  of  the  dislocation  core  geometry  from  electronic 
structure  calculations  is  now  possible  with  current  tools 
and  computational  resources. 
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